Spectrum and transition rates of the XX chain analyzed via Bethe ansatz 
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As part of a study that investigates the dynamics of the s = ^ XX Z model in the planar regime 
|A| < 1, we discuss the singular nature of the Bethe ansatz equations for the case A = {XX model). 
We identify the general structure of the Bethe ansatz solutions for the entire XX spectrum, which 
include states with real and complex magnon momenta. We discuss the relation between the spinon 
or magnon quasiparticles (Bethe ansatz) and the lattice fermions ( Jordan- Wigner representation). 
We present determinantal expressions for transition rates of spin fluctuation operators between Bethe 
wave functions and reduce them to product expressions. We apply the new formulas to two-spinon 
transition rates for chains with up to A'' = 4096 sites. 



I. INTRODUCTION 

The key to a meaningful interpretation of experimen- 
tal or computational data for the low-temperature dy- 
namics of quantum many-body systems is a thorough 
understanding of the nature of the physical vacuum and 
the dynamically relevant collective excitations including 
their quasiparticle composition. In completely integrable 
systems, the quasiparticle configurations that produce 
particular collective excitations can be investigated close 
up. The identity of the former is preserved by conser- 
vation laws notwithstanding their mutual interaction. 
Integrable Hamiltonians that depend on continuous pa- 
rameters make it possible to observe how the physical 
vacuum transforms gradually and, occasionally, changes 
abruptly across a quantum phase transition. Along the 
way, the configurations of quasiparticles in the collective 
excitations are subject to change as well. 

Recently, we investigated the metamorphosis of the 
physical vacuum and the dynamically relevant quasipar- 
ticles (magnons, spinous, ■0; V'*) of tl^c ID s — ^ Heisen- 
berg antiferromagnet {XXX model) by varying the (in- 
tegrability preserving) external magnetic field. ''"^ Intro- 
ducing instead a uniaxial exchange anisotropy also pre- 
serves integrability. The parameter A, which controls the 
anisotropy in the ID s = ^ XX Z model, ' 
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H — ^ {SnS^_^_l -t- S'^S^^i -f AS'^S'^_,_;^} — A— , (1.1) 

n—l 

affects the physical vacuum and the quasiparticle config- 
urations differently. In both situations, the Bethe ansatz 
is an ideal framework for studying quasiparticles, their 
transformations and their interactions. 

The main focus in this paper is on a technical point of 
considerable importance in the study of the XX Z model, 
namely the identification of the general structure of the 
Bethe ansatz solutions of all eigenstates for the case A = 
{XX model). This will facilitate tracking all XX Z 
Bethe ansatz solutions across the planar regime, |A| < 1. 

At A = all states can be characterized as non- 
interacting composites of fermions. Understanding the 



relationship between the magnon, spinon, and lattice 
fermion quasiparticles is important for the interpretation 
of the excitation spectrum via dynamical probes as re- 
alized experimentally or computationally. Furthermore, 
recent advances in calculating transition rates via Bethe 
ansatz'"' offer opportunities to extend the list of 
exact results for dynamical properties of the XX chain. 

In Sec. II we analyze the singularities of the XX Z 
Bethe ansatz equations for A ^ 0. In Sec. Ill we discuss 
the mapping between the spinon and fermion composi- 
tions of the XX spectrum. The two-spinon spectrum of 
the planar XX Z model including the XX limit is dis- 
cussed in Sec. IV. In Sec. V transition rate formulas for 
the XX Z model are introduced and further processed for 
the XX limit. The two-spinon part of the dynamic spin 
structure factor S \^{q,uj) is discussed in Sec. VI. 



II. BETHE ANSATZ EQUATIONS 

We consider Hamiltonian (1.1) for even and periodic 
boundary conditions over the range < A < 1 of the 
anisotropy parameter. In the invariant subspace with z- 
component of the total spin = N/2 — r, all eigenstates 
are represented by r interacting magnons. The magnon 
momenta ki are conserved in the scattering processes and 
are determined by the Bethe ansatz equations:^^ 
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1 _^gi(fe.+fc,) „ 2Ae"''"l 



1 _|_ gi(fci+fcj) _ 2Ae''=J' 



i = 1, 



(2.1) 

The energy and wave number of any such r-magnon state 
depend on the magnon momenta alone: 



E = ^ (cos h~ A) , k = ^h 



(2.2) 



1=1 



At first glance it looks as if Eqs. (2.1) are drastically 
simplified in the limit A ^ with all magnon momenta 
restricted to solutions of e'^'"'* = (—1)' "^. However, a 
closer look reveals that magnon pairs with momenta hi + 
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TT for A 



are a common occurrence. This 



opens the door to nontrivial real and complex solutions 
of Eqs. (2.1). The singular behavior is related to level 
degeneracies. Such level crossings, which also occur at 
other values of A, have been traced back to a realization 
of the sl2 loop algebra symmetry at A = {q+q^^)/2 7^ ±1 
with = 1 for iV > 2.i-i-^'i^ 

In the following, we use the anisotropy parameter 

7 = arccosA (0 < 7 < 7r/2) (2.3) 

and transform the magnon momenta into the rapidities^^ ' 



7 ki 

yi = tan-coty, i = l,...,r. 
The Bcthe ansatz equations (2.1) thus become 



C2yi + i 
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ci(2/i-2/j)+i(l-2/i:2/j) 



- ly V; ciiui - Vj) - i(l - ViV-i) 



(2.4) 



(2.5) 



for i = 1, . . . , r, with ci cot 7 and C2 = cot(7/2). Tak- 
ing the logarithm yields 



Vi ~ Vj 



ViVj 



i = 1 



(2.6) 

with (j){x) = 2arctan(2;). Expressions (2.2) for energy 
and wave number now read 



E = 

C2 ^ + C2 ^ (c22/i) ^ + C2j/j ' 



(2.8) 



The trigonometric Bcthe ansatz equations (2.6) have the 
advantage that each solution is characterized by a set 
of (integer or half-integer) Bethe quantum numbers li. 
These discriminating markers are used to count and clas- 
sify the solutions and are essential for numerical algo- 
rithms designed to find solutions. 



A. XX-limit 

Here we describe the general structure of the solutions 
of the Bethe ansatz equations (2.6) in the limit A — > 0, 
implying ci ^ and C2 1- There exist regular so- 
lutions and singular solutions. The latter are character- 
ized by the occurrence of pairs of rapidities yi,yl with 
the property yiyl = 1. For any such critical pair, the 
argument of </> on the right-hand side of Eqs. (2.6) is in- 
determinate and must, therefore, be treated as a limit 
process. We shall see that some limiting singular solu- 
tions are real while others are complex. An important 
fact is that the simplified structure of the Bethe ansatz 



equations at A = afford a universal treatment of the 
singularities for arbitrary N. 

In the absence of any critical pair of rapidities, 
Eqs. (2.6) decouple and yield the solutions 



Vi = tan 



^ h = n-^I, (2.9) 



for i — 1, . . . ,r, which are all real. The energy of any 



such regular state is 



E^_yIkyi = ycosh 



(2.10) 



B. Real critical pairs 

Now let us assume that among the set of rapidities 
2/1, . . . , there is one critical pair, 

y^jV^, = 1 ^ kj + kj, = TT (mod 27r), (2.11) 

and that it is a real pair. Substituting the ansatz 

y,^y^ + ciS,, z = l,...,r, (2.12) 

into Eqs. (2.6) and taking the limit ci — > then yields the 
non-critical and critical rapidities from successive orders 
in a ci-expansion. The non-critical rapidities yi, i ^ 
arc given by Eqs. (2.9) as in regular solutions. The ci- 
expansion of Eqs. (2.6) leads to the following equation 



(2.7) determining the critical rapidities 2/°, 2/° 



(2.13) 



Here we have introduced 2/^ = 5 (2/? ± 2/?- ) and we use 



^1 ~ AT 



sin kj 



N (2/+)-i -sinfc,' 



(2.14) 



sgn 



N 



(_1)^.-^.*.(2.15) 



The non-critical ki are from (2.9). Criticality implies 
(2/^)^ — [yjY' ~ 1. Note that the contributions of the 
critical rapidities to the energy (2.10) cancel out. Equa- 
tion (2.15) implies that the Bethe quantum numbers of 
the critical pair must satisfy 



(2.16) 



The generalization of this solution to the case of s real 
critical pairs, 



(2.17) 
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with Bethc quantum numbers satisfying I = -^/2 

is straightforward. The non-critical rapidities arc as in 
Eq. (2.9) and the critical pairs are determined from 



1 - 0, 



(2.18) 



with (,ji,o-ji as defined in (2.14) and (2.15). 



C. Complex critical pairs 



D. Generic case 

Suppose we have t critical pairs of complex conjugate 
solutions, yi^ui+ ivi =y^,..., yst-i = + ivf = y2t, 
and s pairs of critical real solutions, J7/ = {jhj*} C {2t + 
1, . . . ,r}, I — 1, . . . , s. Then the real solutions in the limit 
ci — > are of the form 

s 

j/° = tan^/„ ii{l,...,2t}u\Jji, (2.26) 



1=1 



Next we consider the case of a single complex conjugate 
critical pair 



(2.19) 



among the set of rapidities yi, ■ ■ ■ ,yr- Equations (2.6), 
rewritten as real equations, read 



2ci [{yi-u){l- yju) + yjv'^] 
{l-yiu)'^ + {y^vY-c\[{yl-uY + v'^ 

for i = 3, . . . , r and 

2C2M 

1=3 



(2.20) 
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N(p 



2C2V 



1-|-C^(lt2+1'2) 



= 2^(/i+/2)+^[2^/,-A^0(c2y.)], 

1=3 

(2.21) 

4uci(l-u^-v^) 



(l-u2-w2)2 + (2wci)^ 



2ci?;(l-y2) 



3=3 



{l-yju)^ + (yjv)^+cl[{yj~u)^ + v^^ 



(2.22) 



with (p{x) = 2atanh (x) for the critical pair. The XX 
limit is performed along the path u — u'^ — ciu^, v — 
— civ^. The only possible complex solutions that can 
survive the limit ci —^ are critical pairs. The non- 
critical roots arc again as in (2.9). The critical pair of 
rapidities is now determined by the equation 



l + v° 
where 



N 



1-6+1"° 



^ , ut, + v'^ = l, (2.23) 



2 u'^ — sin ki 

" iV ^ (uO)-i -sinfci' 



(2.24) 
(2.25) 



As is custom,^ we shall replace Ii -I- /2 by a single Bethc 
quantum number ± N/2. We then have = for 
all complex critical pairs. 



(2.27) 

with 



r + -70 

j2 ^ * 



For the complex solutions we obtain 
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1 ~ + ^° 

1, i = 1 t 



(2.28) 



(2.29) 
(2.30) 



with 



j=2t+l 



uy — sm k'j 



0\-l 
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Binfc" 



, (2.31) 

Energy and momentum of the state are determined by 
the non-critical roots alone: 



E 



E 

i=2t+l 



cos 



(2.32) 



n{r-s-t)-— /i(mod27r). (2.33) 



N 



i = 2t+\ 



This prescription for handling Bethe ansatz solutions in 
the XX limit will guarantee that all XX Z Bethe eigen- 
states can be traced continuously across the point A = 0. 
Because of the higher symmetry at A = and the asso- 
ciated level degeneracies, the relationship between 
the (real and complex) Bethe eigenstates and the (always 
real) Jordan- Wigner eigenstates of the XX chain is non- 
trivial. The two representations will now be compared 
for all energy levels. 



III. SPINONS VERSUS FERMIONS 

The full spectrum of the XX Z model can be accounted 
for as composites of interacting spinons with spin 1 /2 and 
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semionic exclusion statistics.^'' For the XX case an alter- 
native and equivalent interpretation of the complete spec- 
trum can be established on the basis of noninteracting 
spinless lattice fermions. How are spinon configurations 
related to fermion configurations in single non-degenerate 
eigenstates and in groups of degenerate eigenstates? 

Consider the 2^-dimensional Hilbert space for even 
N divided into subspaces characterized by n+ spinous 
with spin up and n_ spinous with spin down. Equivalent 
quantum numbers are the total number of spinous and 
the z-component of the total spin: 

2n = 7T,_|- + 7i_, 2S^ = 71+ — n_. (3.1) 

The dimensionality of each such subspace is^'' 

da = i(iV + l)-i^K,-<5..0,(3.3) 

a' 

where tr = ± denotes the spinon polarization. Summing 
2n — over n or n+ yields 

N/2 / s 2?i / X 

n=l ^ / ri+=0 ^ ^ 

respectively. The double sum yields 2^. In Table I we 
list the subspace dimensionalities for the case = 8. 



where the allowed values of the fermion momenta pi de- 
pend on whether the number Nf of fermions in the sys- 
tem is even or odd: 

Pi e {{2TT/N){n + 1/2)}, (iV/ even) (3.6a) 
Pi e {{27r/N)n}, {Nf odd), (3.6b) 

for n 0, 1, . . . , — 1. The number of fermions (0 < 
Nf < N) is related to the quantum number in the spin 
representation: = N/2 — Nf. No matter whether Nf 
is even or odd, there are N distinct one-particle states. 
Every one-particle state can be either empty or singly 
occupied, yielding a total of 2^ distinct many-particle 
states with energy and wave number 

Nf Nf 

E - Ef ^'^cospi, k^'^pi. (3.7) 

i=l i=l 

The lowest-energy state in each S'f.-subspace is unique. 
The fermion configuration in reciprocal space of each 
lowest-energy state for iV = 8 is shown in Fig. 1. 
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TABLE I: Dimensionalities W{n+,n-) of the invariant sub- 
spaces with 2n spinons and spin = n+ — n_ for a chain 
of length N — 8. The equivalent quantum numbers in the 
fermion representation are Nf — N/2 — Sf., ric — n — \S^\. 
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The XX Hamiltonian in the fermion representation, 
transformed from (1.1) at A = by the Jordan- Wigner 
mapping to a system of free spinless fermions, reads'' 

i?/ = ^ cos p cj^Cp, (3.5) 
p 



FIG. 1: Configuration in reciprocal space of the Nf — N/2 — 
Sf. fermions in the lowest-energy eigenstate for given of 
a chain with N = 8 spins. The positions of the fermions are 
denoted by full circles, those of vacancies by open circles. The 
fermion momenta pi are in units of 2tt/N. The vertical bars 
are at wave numbers (also in units of 2n/N) — N/i — St /2 
and k+ = 3iV/4 4- Sf./2. 

For the further subdivision of each S'|,-subspace, we 
introduce the wave numbers 

shown as vertical bars in Fig. 1, dividing the band into 
two regions, one in the center, one in the wings. Given 
that the lowest-energy state at S"!, > has all particles 
in the center region, we can characterize all excitations 
as 7ic-particle states. Likewise, all excitations from the 
lowest-energy state at < can be characterized as Tic- 
hole states. The range of this second quantum number is 
ric = 0, 1, . . . , N/2 — \S^\. The number of fermion states 
characterized by ric particles or tic holes in the above 
sense is then found to be expression (3.2) previously used 
for spinons, now with 

7i = 7ic+ IS*!,!, n±^nc + \S^\±S^. (3.9) 
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The next goal is to bring the mapping down to 
energy levels (at fixed wave numbers) within each 
(ri+, n_)-subspace or, in the fcrmion language, (S^^ric)- 
subspace.^^ We note that for given > the num- 
ber of fermions is equal to the number of magnons in 
Bethe ansatz states. Moreover, the fermion energy- 
momentum relation (3.7) is equivalent to the magnon 
energy-momentum relation (2.10). The exception are the 
critical magnon pairs, whose momenta are different (real 
or complex) from the corresponding fermion momenta 
(always real). Critical pairs only occur in degenerate 
levels and do not contribute to the energy of the state. 

This mapping between energy levels provides a use- 
ful tool for determining the Bethe quantum numbers of 
eigenstates directly from the momenta of the fermion 
configurations. For all non-critical rapidities, we have 

^h=n-p,. (3.10) 

For real critical pairs, the Bethe quantum numbers are, 
in general, shifted relative to the positions predicted by 
(3.10), but in such a way that \Ij + Ij* \ = N/2 is main- 
tained. Complex critical pairs are specified by a single 
Bethe quantum number /^*^ = as discussed in Sec. II C. 
In the following, we look more closely at the two-spinon 
excitations in the two representations. 

IV. TWO-SPINON SPECTRUM 
A. Planar XXZ model 

Returning to the XXZ model (1.1), we note that there 
exist two-spinon states with = 0, ±1. At A = 1 
these states are either triplets (with total spin 5t = 1) 
or sing lets {St = 0). There are lN{N + 2) triplet 
levels and IN{N - 2) singlet levels, distinguishable by 
their Bethe quantum numbers.^"'"" Integrability guaran- 
tees that each eigenstate anywhere in the planar regime 
can be traced back to the point of higher symmetry with- 
out ambiguity. This justifies that we use the St multiplet 
labels for states at A < 1. 

At finite N and A < 1, the two-spinon triplet compo- 
nents with S^ = ztl remain degenerate (because of reflec- 
tion symmetry in spin space) but are no longer degener- 
ate with the two-spinon triplet components with S^ = 0. 
We shaU see that at A = 0, a subset of ^N{N - 2) 
states from the latter set become degenerate with the 
two-spinon singlets. 

We have tracked the Bethe ansatz solutions of every 
two-spinon state for = 8 from A = 1 to A = 0, in- 
deed across this highly singular point of the Bethe ansatz 
equations all the way to A = —1. Here we briefly report 
on what we found along the stretch < A < 1.-^ For 
all states, solutions of the Bethe ansatz equations were 
found for which all magnon momenta ki (or all rapidities 
Hi) vary continuously across the parameter regime. 



The two-spinon triplets with S^ = 1 are easy to handle 
numerically. All rapidities are finite and real. Two- 
spinon triplets with S^ = 0, by contrast, tend to pose 
some computational challenges. At A = 1 all such states 
start out with one fc^ = magnon, implying yi = ±oo. 
While this singular behavior is benign,-" numerical prob- 
lems are caused by the fact that in some states yi stays 
infinite at A < 1, whereas in other states it becomes fi- 
nite. The majority of states from this set have one real 
critical pair at A = 0. To ensure continuity of the Bethe 
ansatz solutions when A is varied, it is sometimes neces- 
sary to change the values of two Bethe quantum numbers 
leaving their sum invariant. 

The two-spinon singlets are far more difficult to handle 
numerically. Each such state has one complex-conjugate 
pair of rapidities. A numerical analysis of two-spinon sin- 
glets at A = 1 was reported in Ref. 20. Again, continuity 
of the solutions may make it necessary to change some 
Bethe quantum numbers along the way. At A = each 
state from this set with one complex-conjugate critical 
pair becomes degenerate with a two-spinon triplet state 
that has one real critical pair. 



B. XX limit 

According to the scheme developed in Sec. Ill, the two- 
spinon states with S^ = 1 comprise all configurations in 
which no state with positive energy in the fermion band 
is occupied. These are the fermion configurations with 
ric = 0. The N/2 + 1 single-particle states which are 
accessible under this restriction and among which the 
N/2 — 1 fermions can be distributed does indeed produce 
N{N + 2)/8 distinct eigenstates. All possible configu- 
rations for ~ 8, generated systematically across the 
allowed range in reciprocal space, are depicted in Fig. 2. 
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FIG. 2: Set of N{N + 2)/8 two-spinon states with St = 1 at 
A = for TV = 8 as specified by the fermion momenta (top 
scale) or Bethe quantum numbers (bottom scale). The wave 
numbers k and fermion momenta pi are in units of 2'k/N. 
The dashed lines represent and fcj. The energies are e = 
E — Eg relative to the ground state with Eg = —2.613. 
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The magnon momenta, rapidities, Bcthe quantum 
numbers, and energies of these states, now sorted accord- 
ing to wave numbers, are Hsted in Tabic II. No critical 
pairs can be formed in any of these states, which makes 
it straightforward to calculate transition rates for fairly 
long chains (see Sec. VI). 

TABLE II: Wave number, energy, magnon momenta, Bethe 
quantum numbers, and rapidities of the A^(A'^ + 2)/8 two- 
spinon triplet components with = 1 at A = for N = 8. 
The quantities k,ki are in units of 2-k/N. The ground-state 
energy is Eg ~ —2.613. 
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The two-spinon states with = comprise all {ric = 
l)-states relative to the ground state configuration in 
the fermion band. Removing one of N/2 fermions from 
the ground-state configuration and placing it into one of 
the N/2 empty single-particle states yields N^/A distinct 
eigenstates but only N{N + 2)/8 distinct levels in the 
(A:, i?)-plane because of degeneracies. All possible con- 
figurations for iV = 8 including the ground state, gener- 
ated systematically across the allowed range in reciprocal 
space, are depicted as rows of full circles in Fig. 3. 

There are N{N — 2)/8 twofold degenerate levels in the 
(fc, i?)-plane with the two states distinguished by one pair 
of fermions. The two fermions in question (full circles 
connected by a line in Fig. 3) add up to tt mod(27r), 
thus contributing nothing to the energy (3.7). Addition- 
ally, there exist N/2 non-degenerate (ric = l)-statcs at 
wave- numbers Nk/2-K = 1, 3, . . . , A'^ — 1 representing the 
highest-energy two-spinon state for those wave numbers. 

In summary, the N"^ / A {ric ~ l)-states in the fermion 
representation at A = originate from the set of 
A^(A^ -I- 2)/8 two-spinon triplets with = Q and the set 
of A^( A^ — 2) /8 two-spinon singlets. Each singlet becomes 
degenerate with a triplet at A = 0, while A^/2 triplets re- 
main non-degenerate. When analyzed via Bethe ansatz, 
the two states of each degenerate two-spinon level with 
= have one critical pair of magnon momenta as 
described in Sec. II. Our analysis shows that the criti- 
cal pair is real for the triplet state and complex for the 
singlet state. This is consistent with the known 2-string 



nature of the two-spinon singlets and 1-string nature of 
the two-spinon triplets at A = 1. 

In the states with real critical pairs, a solution may 
only exist for Bethe quantum numbers that are shifted 
relative to those predicted by (3.10) but still satisfy 
(2.16). The latter are indicated in Fig. 3 by open cir- 
cles. In states with complex conjugate critical pairs, two 
of the Bcthe quantum numbers (3.10) arc replaced by a 
single number I^*^ = representing the critical pair. 

The actual Bethe ansatz solutions for A^ = 8 pertain- 
ing to all A^^/4 two-spinon states with = Q are listed 
in Table III (triplets) and Table IV (singlets). In some in- 
stances, different configurations of Bethe quantum num- 
bers lead to equivalent Bethe wave functions. Only one 
configuration is represented in Fig. 3.^^ 
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4 1 • * * ! * o £ = 0.765 

2 I ' * ' O^r^ * £=l-307 

3 ; * * * ; * £=1.848 

4 1 • g * j * £=1.848 

5 o S " • £ = 1.307 



1; 4 3 2 1 -1 -2 -3 -4 

FIG. 3: Ground state and set of N'^/A = N(N + 2)/8 + 
N{N - 2)/8 two-spinon states with = for A = 8. The 
positions of the full circles in relation to the top scale denote 
the fermion momenta and those in relation to the bottom scale 
the Bethe quantum numbers predicted via Eq. (3.10). The 
critical pairs among them are connected by a solid line. The 
open circles and squares denote the actual Bethe quantum 
numbers associated with the critical pairs for real and complex 
rapidities, respectively. The wave numbers k and fermion 
momenta pi are in units of 2tt/N. The dashed lines represent 
kc and k^ . 
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TABLE III: Wave number, energy, magnon momenta, Bethe quantum numbers, and rapidities of the A''(A'^ + 2)/8 two-spinon 
triplet components with = at A = for N = 8. The quantities k, ki are in units of 2-k/N . 



K t, — rjG 


ki k2 ks ki 


T T T T 

11 12 13 14, 


yi y2 yz yt 


0.000 


2.5 3.5 4.5 5.5 


1.5 0.5 -0.5 -1.5 


0.668 0.199 -0.199 -0.668 


1 0.765 


2.5 3.5 4.5 6.5 


1.5 0.5 -0.5 -2.5 


0.668 0.199 -0.199 -1.497 


2 1.307 


2.5 3.5 5.110 6.890 


1.5 0.5 -1.5 -2.5 


0.668 0.199 -0.466 -2.147 


3 1.307 
1.848 


2.5 4.5 4.728 7.272 
2.5 3.5 5.5 7.5 


1.5 -0.5 -0.5 -3.5 
1.5 0.5 -1.5 -3.5 


0.668 -0.199 -0.294 -3.402 
0.668 0.199 -0.668 -5.027 


4 0.765 
1.848 


3.5 4 4.5 
2.5 4 5.5 


0.5 -0.5 -0.5 -3.5 
1.5 -0.5 -1.5 -3.5 


0.199 -0.199 -oo 
0.668 -0.668 -oo 


5 1.307 
1.848 


0.728 3.272 3.5 5.5 
0.5 2.5 4.5 5.5 


3.5 0.5 0.5 -1.5 
3.5 1.5 -0.5 -1.5 


3.402 0.294 0.199 -0.668 
5.027 0.668 -0.199 -0.668 


6 1.307 


1.110 2.890 4.5 5.5 


2.5 1.5 -0.5 -1.5 


2.147 0.466 -0.199 -0.668 


7 0.765 


1.5 3.5 4.5 5.5 


2.5 0.5 -0.5 -1.5 


1.497 0.199 -0.199 -0.668 



TABLE IV: Wave number, energy, magnon momenta, Bethe quantum numbers, and rapidities of the A''(A'^ — 2)/8 two-spinon 
singlets at A = for N = 8. The quantities k, ki are in units of 2tt/N. 



k E -Eg 


k* k:i ki 


h h 


u V yz yi 


2 1.307 


4 3.5 2.5 


0.5 1.5 


-0.765 -0.644 0.199 0.668 


3 1.307 


4 4.5 2.5 


-0.5 1.5 


-0.541 -0.841 -0.199 0.668 


4 0.765 
1.848 


4 4.5 3.5 
4 5.5 2.5 


-0.5 0.5 
-1.5 1.5 


1 -0.199 0.199 
1 -0.668 0.668 


5 1.307 


4 5.5 3.5 


-1.5 0.5 


0.541 0.841 -0.668 0.199 


6 1.307 


4 5.5 4.5 


-1.5 -0.5 


0.765 0.644 -0.668 -0.199 



V. MATRIX ELEMENTS VIA BETHE ANSATZ 
FOR A = 



We start from the determinantal expressions for the ^'^>' ^'^^ l^r{{yi}) J det K({j/i}) det K({y°}) ' 



ICr■±li{y^})\^' |dctr±| 



transition rates 



where 



= E M = ^.+^-- (5-2) 



X G(jj'^ yt) 

These expressions, which were derived in Ref. 10 based ^ab = FN{y^Tyb) | „, q r JJ^ ' 

on previous work reported in Refs. G,7,8 have the form 



G{yly,) \}^G{y,,y,) 



detK({yn)detK({,J)' G*(,o, 11 y,) 



(5.4) 



between XX Z eigenstates characterized by real ki for '^ab — \ cos7'^ K{ya,yj) ^ ' 

the spin operators ^ ]y~ Kijja) ■ o, — b 



(5.6) 
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ab 



I n(r,, ^ nin, ,,0^ Switching back from the non-critical rapidities via (2.9), 

^N(ya,yb) I ~7y o\~ 11 r'( o — oT ~ cot(/Ci/2), to the non-critical magnon momenta ki, 



we can bring expression (5.4) into the form 



G*(,.,,0) n I ' (5-7) n n cos ^ 



Cr+i = 1> a= l,...r+l, 6= 1, 



M+{q) = —^zll^ |dct5+|2, (5.19) 

n C0S2 ^ n C0S2 ^ 
i<j i<j 



■ — 1 . . . r 



'^i, = S sin /cq - sin fcP , a = 1, . . . , r + 1. 



G (ya. yb) fJi G*{yj,yb) J ^ 
r- = 1, a= l,...r, 6= l,...,r-l. n cos^ ^ cos^ 



M^{q) = — i |det5-r, (5.20) 



2 



^r({yj) = ni^(y-2^^)l' 



2/N 



I • 5 = 1 r — 1 

= < sin fcO - sin fcfc ' ' ' ' , a = 1, . . . , r, 



1 



K( (l-i/^)(l-?/-^)sin^ 

K[y,y) = z ^tt; „x . 2 ■ (^-lO) 



i^iy) 



(y-y')' + (l-y')(l-y")sin'7 



(y ^ - ;/)sin7 
2;cot(7/2) + (ycot(7/2))- 



G(y,y') 



{y-y')cotj + i{l-yy') 



(5.11) 



(5.12) 



where are Cauchy-type matrices, whose determinants 
can be evaluated explicitly: 



r+l 



det5+ 



n (sinfcf - sinfc°) JJ (sinfcj - sinfe;) 

2 \ i<j i<j 



N 



dct5" 



N 



n'n(sinfc,-sinfcO) 

j=ii=i 

r r—1 



n'n(sinA:0-sinA:,) 



i^jv(y,2/') = 



1 v^i - 1 + + (?/2 - 1) C0S7 



27V (y~y')sm7 

(5.13) 

Performing the XX limit, 7 — > 7r/2, in these expressions 
is straightforward as long as no critical pairs of rapidities 
are present: 



K{y.y') 



i^iy) 



jl-y^Kl-y'^) 
(1 - yy'y 

y^^ - y 
y +y 



G{y,y') = i 



1 - yy' 



FN{y,y' 



v/i-yVi-2/" 
1 yT^?2i + y2 



(5.14) 
(5.15) 
(5.16) 

(5.17) 



2N yr^^ y-y'' 
detK = 1. (5.18) 



This reduces the transition rates for the perpendicular 
spin fluctuations (between states without critical pairs) 
to compact product expressions: 



M+{q) = 



r+l 

n 



n sm n sm 



r r+l „ 

UN^ U sin' 

77'' 1.0 ^o 

n sin' n sin' ^ 
Ma^(9) = 

'■-1 1,0 7 

n^v'nsin'^ 

i=i i=l 



(5.21) 



(5.22) 



In the corresponding reduction of the transition rate 
(5.3) for the parallel spin fluctuations, a complication 
arises, caused by the possibility that some magnon mo- 
menta of the two states might be identical. However, 
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this singular behavior turns out to be instrumental for 
the exact evaluation of M^{q). A nonzero result is only 

possible if the two sets of Bethe quantum numbers {/p^} 
and {li} differ by no more than one clement. For all such 
transitions the rate is 

M|(<?) = ^ (5.23) 

in agreement with a well-known result derived in the 
fermion representation.^ ' 

In the following application of the transition rate ex- 
pressions to a T = spin dynamic structure factor of 
the XX model we have chosen a situation where excited 
states without critical rapidities arc important. The cal- 
culation of transition rates with critical pairs requires 
further developmental work. 

VI. TWO-SPINON TRANSITION RATES 

From recent studies in the framework of the algebraic 
analysis for the infinite chain,-"''-"' we know that the rel- 
ative integrated intensity of the two-spinon contribution 
to the dynamic structure factor 

S.+ {q,Lo)^2TrY,^^xiQ)H^-^x) (6-1) 

A 

probing the spin fluctuations perpendicular to the sym- 
metry axis of the XX Z model is 73% for the Heisenberg 
case (A = 1) and steadily growing toward 100% on ap- 
proach of the Ising case (A = oo). 

The nonzero two-spinon intensity is reflected in the 
reciprocal finite- TV scaling behavior of the transition rates 

(2) — 

M_^{q,uj\) = NM-y (q) and the scaled density of states 

D^^\q,iOx) - 2V[iV(wA+i - ujx)], (6.2) 
which makes the product 

S^^l{q,uj) = M^_'l{q,uj)D^^\q,uj) (6.3) 

converge toward a picccwisc smooth function in the limit 
iV — > (X), representing the result of the infinite chain. 

A qualitatively different finite-A'^ scaling behavior is 
found in the XX case (A = 0) for the two-spinon tran- 
sition rates contributing to S \.{q,uj). For the tran- 
sition rates to converge toward a non-vanishing piece- 
wise smooth function they must be scaled differently: 
M^^l{q,oj) = N^/'^M-{q). This is illustrated in Fig. 4 
for q = n. In the main plot we show data for N ~ 
512, 1024, 2048, 4096 of iV3/2M^(7r) versus wa. The scal- 
ing is near perfect across the band. The divergence build- 
ing up as TV ^ oo in this quantity is stronger, ~ i^^'^, 
than the known infrared singularity in the dynamic struc- 
ture factor,-'' S ^{q,Lo) ~ as documented by the 

inset to Fig. 4. 



Given the non-reciprocal scaling behavior of the tran- 
sition rates and density of states, the relative intensity of 




1 2 

<»>. 

FIG. 4: Scaled two-spinon transition rates M^}^{tv,uj) = 
A^/^Af^(7r) at A = (data for chains of size N = 
512,1024,2048,4096). The inset shows the same quantity 
multiplied by ujf. 



the two-spinon dynamic structure factor S_!^{q,uj) van- 
ishes in the limit A^ — > oo. Hence the singularity struc- 
ture of the hmction M_^{q,uj) has no direct bearing on 
the singularity structure of S \.{q,uj). A distinct sin- 
gularity structure and spectral-weight distribution which 
is a property of all 2m-spinon excitations combined will 
emerge in the limit N —f oo. 

Consequently, the exactly known leading singularities 
at w = 0, 1, 2, ... in the frequency-dependent spin auto- 
correlation function $Q+(tj) = J^^ {dq/2Tr)S ^{q^uj), as 

worked out in Ref. 27, for example, arc not attributable 
to specific 2m-fermion excitations, because the integrated 
intensity of each 2TO-fermion contribution taken in isola- 
tion is likely to vanish in the limit N —^ oo. 

To reconstruct the spectral-weight distribution of 

S |_((7, w) at A = and to determine its singularity 

structure in the limit A^ ^ oo from finite-TV data for 
excitation energies and transition rates we need to be 
able to properly handle Bethe wave functions with criti- 
cal pairs of rapidities. We already know (Sec. II) how to 
solve the Bethe ansatz equations for all eigenstates in the 
limit A — > oo. One challenging problem for the calcula- 
tion of transition rates is that Bethe wave functions with 
critical pairs thus obtained vanish identically as pointed 
out in Ref. 2<s. However, our numerical analysis strongly 
suggests that the vanishing norm HV'aII in the denomina- 
tor of (5.1) is compensated by the vanishing transition 
rate in the numerator to produce a unique finite ratio in 
the hmit A ^ 0.-'^ 
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